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Semileptonic decays of D mesons in unquenched lattice QCD* 
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We present our preliminary results for semileptonic form factors of D mesons in unquenched lattice QCD. 
Simulations are carried out with rif =2 + 1 dynamical quarks using gauge configurations generated by the MILC 
collaboration. For the valence quarks, we adopt an improved staggered light quark action and the clover heavy 
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1. INTRODUCTION 

In order to extract the CKM matrix element 
|Kt&| accurately with experimental measurement 
of the semileptonic decay width, a precision lat- 
tice QCD calculation of the B — > ir form fac- 
tor is required. To check the reliability of lattice 
calculations of the heavy to light form factors, 
study of the semileptonic decays of D mesons, 
such as D — > K and D — ► tt, is a good test 
ground because the corresponding CKM matrices 
\V CS \ and \V c d\ are known more accurately than 
\V u b\. Furthermore, forthcoming experiments by 
the CLEO-c collaboration will provide more strin- 
gent checks of lattice calculations in the D meson 
system. 

We have started new lattice calculations of 
heavy quark physics in unquenched QCD F, and 
here we report our preliminary results for semilep- 
tonic form factors of D mesons. We are us- 
ing unquenched gauge configurations with rif = 
2 + 1 improved staggered quarks generated by the 
MILC collaboration^, with which the systematic 
errors due to the quenched approximation should 
be almost absent. For the valence light quarks, 
we adopt an improved staggered quark action, 



which allow us to simulate at lighter quark mass 
than previous studies with the Wilson-type light 
quarks. Hence, our new calculations should have 
a better control over the chiral extrapolations. 

2. METHOD 

In order to combine the staggered light quark 
with the Wilson-type heavy quark in heavy-light 
bilinears, we convert the staggered quark prop- 
agator g(x, y) to the "naive" quark propagator 
G(x, y) according to 

g(x,y) n(x) nt(„) = G(x,y) (1) 

with Ofx) = 7o°7r72 2 73 3 |Sl' The 3-point func- 
tion for the matrix element is then computed as 



C 3,7* r (*»»*y;Pn->P£») 



o*(PD-P7r)y-ipDX 



(Tr[gl(y,0)n\y) l5%l G c (y,x) l5 n(x)g u (x,0)]), 

where subscripts d, c, u denote quark flavors. The 
matrix element can be extracted from the ratio 
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with the D meson (Wilson-Naive) 2-point func- 
tion and the pion (Naive-Naive) 2-point func- 
tion CJ. Care is needed, however, for the overall 
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Table 1 

Quark mass, statistics and the sink time. 
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normalization of amplitude since the naive quark 
action describes 16 fermions, which can cause the 
doubling of these correlation functions. 

In Ref. |3] it is shown that the Wilson-Naive 
2-point function does not have the doubling 
because contributions of quarks with momentum 
p ~ 0(n/a) are suppressed by the Wilson term. 
The same also holds for the 3-point functions 
which include at least one Wilson propagator such 
as . On the other hand, the Naive-Naive 2- 

point function CJ should have 16 equivalent con- 
tributions. Therefore one has to divide it by 16 
to get the physical amplitude; ,phys = CJ/16. 

3. SIMULATION 

Unqucnchcd calculations are performed using 
rif = 2 + 1 dynamical gauge configurations ob- 
tained with an improved staggered quark action 
on a 20 3 x 64 lattice (a' 1 w 1.58 GeV) 2 . For 
the valence light quarks we use the same stag- 
gered quark action as for the dynamical quarks. 
The valence light quark (u, d) mass m ral is usu- 



ally set equal to the dynamical light quark mass 
TO soa p Qr va i ence c harm quark we use the 
clover action with the Fermilab interpretation^. 
The hopping parameter is fixed to -fQharm = 
0.119, based on our spectrum study [I]. The 0(a) 
rotation^] is performed for the vector current. 

The 3-point functions are computed in the D 
meson rest frame (pn — 0) for the light meson 
momentum up to (1, 1, 1) in lattice units, us- 
ing local source and sink. The sink time is fixed 
to t x = 20 — 26 depending on to™ 1 , whereas the 
source time is set to to — with an exception at 
to™ 1 = 0.01, where we average over results from 
four source times to = 0, 16, 32 and 48. Some sim- 
ulation parameters are summarized in Tabled 

In addition to the unquenched calculations, we 
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Figure 1. D s — > r) s (ss) form factors. 

also perform a quenched simulation at To; = 
to™ 1 = 0.0415 using /3 = 5.9 (a^OO GeV) con- 
figurations on a 16 3 x 32 lattice used in our previ- 
ous study Comparison between the quenched 
result with the staggered light quarks and that 
with the Wilson- type light quarks |S| allows us to 
check the validity of our new calculations. 

For the vector current renormalization Zy 1 



we follow the method in Ref. We take 

Z% = PVik (ZfZffl\ where Z q v q (q = c,d) 
is the renormalization constant for the flavor- 
conserving current, which we compute nonpertur- 
batively from the charge normalization condition 
Z v q (D(0)\V^ q \D(0)) = 2m D . The p Vfi is set to 
unity. The one-loop calculation is in progress. 

4. RESULTS 

Form factors are defined through 



MV*\D) = f + (q 2 



PD + Prr 



+ fo(q 2 ) ■ 



= V2^ [v»f {] (E)+ptf ± (E)] 

with q — p D -pn, v = pD/m D , p± = p^-Ev and 
E = En. The second expression using /|| and f± 
is more convenient when one considers the heavy 
quark expansion and the chiral limit. 

4.1. D s ^r) s 

In Fig. ^ we summarize the results of form 
factors fo and /+ for the D s — > rj s (ss) decay 
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Figure 2. Chiral extrapolation for f\ 
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Figure 3. D — + tt and D K form factors. 

obtained with the naive (staggered) light quarks 
and previous results with the Wilson- type light 
quarks. The lines in the figure are fitting curves 
with a parametrization by Becirevic and Kaidalov 
(BK) 6 . One can see a nice agreement be- 
tween the quenched result with the naive quarks 
(squares, dashed line) and that with the Wilson- 
type quarks (circles, dotted), showing that our 
new method works well. We also note that the un- 
quenched result (diamonds, solid) is larger than 
quenched ones for q 2 > 0. See also Ref. J7J for a 
similar comparison for the B s — > i] s form factors. 

4.2. D -> tt and D -> K 

To obtain D — > ir/K form factors at the phys- 
ical quark mass, we need to perform a chiral ex- 
trapolation using data in range of to™ 1 = 0.01- 
0.03. We do this for /|| and f± at fixed pion(kaon) 
energies E„ik) because the chiral perturbation 



formulas for the heavy to light form factors are 
given in such a wayjHj- In order to interpolate 
and extrapolate the results to common values of 
we use a fit with the BK parametrization. We 
then perform a linear chiral extrapolation in to™ 1 
at nine values of (aE^) 2 . One example of these 
procedures is shown in Fig.|21for f®~ >7r . Finally 
/|| and /j_ are converted to fo and /_)_. 

The D — > 7r and D — > K form factors are shown 
in Fig. |3| together with experimental values at 
q 2 =00. Our results at q 2 = are 



f^ K (0) = 0.75(3), f^'(0) = 0.64(3) (3) 

with statistical errors only, whereas experimental 
values are f^ K (0) = 0.73(2) and f^ w (0) = 
0.73(13) with \V CS \ = 0.996(13) and \V cd \ = 

0. 224.16) 9^. Our results are in agreement with 
the experimental values. The analysis including 
the chiral logarithm and the one-loop renormal- 
ization constant is underway. 
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